Wave modes trapped in rotating nonlinear potentials 
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We study modes trapped in a rotating ring with the local strength of the nonlinearity modulated 
as cos (20), where is the azimuthal angle. This modulation pattern may be of three different 
types: self- focusing (SF), self-defocusing (SDF), and alternating SF-SDF. The model, based on 
the nonlinear Schrodinger (NLS) equation with periodic boundary conditions, applies to the light 
propagation in a twisted pipe waveguide, and to a Bose-Einstein condensate (BEC) loaded into a 
toroidal trap, under the action of the rotating nonlinear pseudopotential induced by means of the 
Feshbach resonance in an inhomogeneous external field. In the SF, SDF, and alternating regimes, 
four, three, and five different types of stable trapped modes are identified, respectively: even, odd, 
second-harmonic (2H), symmetry-breaking, and 2H-breaking ones. The shapes and stability of these 
modes, together with transitions between them, are investigated in the first rotational Brillouin zone. 
Ground-state modes are identified in each regime. Boundaries between symmetric and asymmetric 
| modes are also obtained in an analytical form, by of a two-mode approximation. 

PACS numbers: 42.65.Tg; 03.75.Lm; 47.20.Ky; 05.45.Yv 

I. INTRODUCTION 

Optical and matter waves exhibit a plenty of dynamical scenarios under the action of effective nonlinear potentials 
(which may sometimes be combined with usual linear potentials) The dynamics of such systems is governed by 
the nonlinear Schrodinger equation (NLSE) in optical media, or Gross-Pitaevskii equation (GPE) in the context of 
Bose-Einstein condensates (BECs). In either case, the nonlinear pseudopotential [2j may be induced by a regular 
0, H| or singular spatial modulation of the local nonlinearity. These systems have been studied theoretically in a 
variety of one- [JIBI and two-dimensional (ID and 2D) Q settings, and recently reviewed in Ref. [lj. To the same 
general class belong models which predict that, in any dimension D, stable fundamental and vortex solitons can be 
supported by a purely self-defocusing (SDF) nonlinearity growing towards the periphery (r — >• oo) at any rate faster 
than r D @. 

£C) ' In optics, such nonlinear potentials may be designed using the mismatch between the nonlinearity of the host 
material and solid [7] or liquid [s| stuff filling voids of photonic-crystal-fiber waveguides. Another possibility to 
create the effective nonlinear potential in optics is offered by inhomogeneous distributions of dopants which induce 
\ the resonantly enhanced nonlinearity [9|. In particular, it is possible to use the Rhodamine B dopant added to the 
SU-8 polymer (a commonly used transparent negative photoresist) fl(jj |. or Pr 3+ ions doping the Y2SiOs host medium 

In BEC, the pseudopotential can be created with the help of the Feshbach resonance controlled by nonuniform 
1 magnetic [l3| or optical Q fields. In particular, the necessary pattern of the spatial modulation of the scattering 
length, which determines the local strength of the cubic nonlinearity in the respective GPE, can be induced by 
appropriately designed magnetic lattices jl5j ]. 

Another well-known tool for the creation of various dynamical states is provided by rotating potentials, which may 
trap optical and matter waves. Effects of the rotation have drawn a great deal of attention in the studies of BEC. A 
well-known results is that rotational stirring of the condensate with repulsive interactions leads to the formation of 
vortex lattices, see review [16|. Under special conditions (the compensation of the trapping by the centrifugal force), 
giant vortices can be produced too, which were studied in detail experimentally [ItJ and theoretically [l8j|. In a binary 
immiscible BEC, vortex streets were predicted to form, instead of the vortex lattices [l9|. On the other hand, it 
was predicted that the rotation of self-attractive condensates gives rise to several species of localized modes, such as 
solitary vortices, mixed- vorticity states ("crescents"), and quasi-solitons [20j . 

It is natural too to consider the dynamics of matter- wave modes trapped in rotating lattices, which can be created 
in BEC by broad laser beams transmitted through a revolving mask [21[. Quantum states and vortex lattices have 
been studied in this setting 12211 , as well as the depinning of trapped solitons and solitary vortices when the rotation 
rate exceeds a critical value [23j. The nucleation of vortices in the rotating lattice was demonstrated experimentally 

In optics, a setting similar to the rotating lattice can be realized in twisted photonic-crystal fibers [25|, [26|. In 
plain optical fibers, the twist affects the polarization dynamics [27j and couples it to the transmission of temporal 
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FIG. 1: (Color online) The pipe nonlinear waveguide with the intrinsic nonlinear potential, twisted at rate uj 



solitons [29|| . In helical photonic -cry stal fibers, modified Bragg reflection and enhancement of the mode conversion 
and transport have been studied [28| , and the transformation of the linear moment of photons into the orbital angular 
momentum has been demonstrated recently [26j |. 

The simplest version of the rotating lattice is represented by the revolving quasi-lD double- well potential (DWP). 
It gives rise to azimuthal Bloch bands [30], and allows one to support solitons and solitary vortices even in the case 
of the SDF nonlinearity [3l|. The generation of a vortex lattice in the rotating DWP was studied too [32]. Further, 
it is well known that the interplay of the DWP with the SF or SDF nonlinearity provides for the simplest setting for 
the study of the spontaneous symmetry breaking of even and odd states, respectively, in one dimension [33|. In this 
connection, a natural problem, which was recently considered in Ref. [UJ], is a modification of the symmetry-breaking 
scenarios in a rotating ring carrying the DWP potential, along with the nonlinearity. 

While the dynamics of nonlinear waves trapped in rotating linear potentials has been studied in detail, previous 
works did not tackle modes pulled by rotating nonlinear (pseudo-) potentials. This setting, which may be imple- 
mented in optics and BEC alike, is the subject of the present work. To analyze the basic features of the respective 
phenomenology, we here concentrate on the ID nonlinear potential on a rotating ring, as shown in Fig. [TJ In optics, 
this system is realized as a pipe (hollow) waveguide with an azimuthal modulation of the local nonlinearity, a (6) (0 
is the angular coordinate), twisted with pitch 2tt/uj, where uj plays the role of the effective rotation speed. In BEC, 
a similar setting corresponds to a toroidal trap, which is available in the experiment [35] , combined with the rotating 
nonlinear potential, that can be superimposed on top of the trap [36|. This combination realizes a rotating ring [37] 
carrying the nonlinear potential. 

It is necessary to define the sign of the nonlinear potential. In this work, we consider three distinct cases, namely, 
the SF (attractive), SDF (repulsive), and SF-SDF (alternating) nonlinearities, all subject to the harmonic spatial 
modulation, which is defined as follows: 

<tsf(0) = - sin 2 (9, (la) 
<tsdf(0) =cos 2 #, (lb) 

^SF-SDF (0) = COS (20). (lc) 

The solution domain is set as — tt < < +7r. Note that, in all the cases the local nonlinearity coefficient (pQ) has its 
maxima at points = and = ±7r, and minima at = ±7r/2. We aim to find basic types of modes trapped in 
rotating nonlinearity profiles (pQ), and establish their basic properties, such as symmetry /asymmetry, stability, and 
the identification of the respective ground states, varying two control parameters, viz., the rotation speed, uj, and the 
total power (norm), P, which is defined below in Eq. (|20|) . 

The paper is structured as follows. In Section II, we formulate the system and methods of the stability analysis of 
the modes. In Section III, we present numerical results for the basic modes and their stability in each type of the 
nonlinear potential (pQ). In Section IV, we present analytical results, obtained by means of a two- mode approximation, 
which explain a boundary between the symmetric and asymmetric modes in each case. The paper is concluded by 
Section V. 
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II. THE MODEL 



The dynamics of the optical wave (or the BEC wave function) in the rotating ring is governed by the normalized 
one-dimensional NLSE (GPE), subject to the periodic boundary conditions: 
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-0(0) = i/j(6 + 27r), where z is the propagation distance in the case of the optical waveguide, and the radius of the ring 
is scaled to be 1. It is more convenient to rewrite Eq. ((2]) in the rotating reference frame, with 9' = — uoz: 
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while the solution domain is denned as above, i.e., — n < 9' < +n. For the BEC trapped in the rotating potential, 
the respective GPE differs by replacing z with time t. Equation (j3J) conserves the total power (norm) of the field and 
its Hamiltonian (energy), 

r + 7T 
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with the asterisk standing for the complex conjugate. Stationary modes with real propagation constant — /i (in terms 
of BEC, \i is the chemical potential) are sought for as 

iP(0',z) = exp(-i f iz)<P(e'), (6) 

with complex function <\> (6 f ) obeying equation 



d6' 



(7) 



The periodic boundary conditions, a(6' + 2tt) — (j{0') and ^{0' + 2tt) = ip(Q'), make Eq. (j3j) invariant with respect 
to the boost transformation, which allows one to change the rotation speed from uj to uj — N with arbitrary integer N: 



(#', z\ uj - N) = (<9', z; uj) exp 



-iNO' 



Nuj\z 



(8) 



hence the speed may be restricted to interval < uj < 1. Furthermore, Eq. (j3]) admits an additional invar iance, 
relating solutions with opposite signs of the speed: i(j(6 f ,z;uj) = ^(0\ —z; —uj). If combined with boost uj —> uj + 1, 
the latter transformation demonstrates that the solutions with 

uj = 1/2 ± 5, (9) 

where S < 1/2, are tantamount to each other. Thus, the rotation speed may be eventually restricted to the fundamental 
interval, 



< uj < 1/2, 



(10) 



which plays the role of the first rotational Brillouin zone, cf. Ref. [30| . 

The dynamical stability of the stationary solutions has been investigated via numerical computation of eigenvalues 
for small perturbations, and verified by direct simulations of the perturbed evolution. The perturbed solutions are 
introduced as usual, 



iP = e-^ z [(t){e') +eu{0')e iXz +^*(<9> 



(ii) 



where e is an infinitesimal amplitude of the disturbance, u{0') and v(6 r ) are the corresponding eigenmodes, and A the 
eigenfrequency. The substitution of ansatz ([TT)) into Eq. ([3]) and linearization leads to the linear eigenvalue problem, 
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where h = —(l/2)d 2 /d(0') 2 + 2a\(j)\ 2 is the single-particle Hamiltonian. The underlying solution <j) is stable if Eq. 
([T2)) generates solely real eigenvalues. Lastly, stationary equation ([7j) was solved using numerical code "PCSOM" 
borrowed from Ref. I38I. 



TABLE I: Different species of stationary modes, labeled by input waveforms which generate them as solutions of Eq. ©. 



Inputs 




Types of modes 


cos<9' 




Symmetric (even) 


sin<9' 




Antisymmetric (odd) 


sin 2 0' 




Second-harmonic (2H) 


b + cos i 


V, < b < 1 


Symmetry-breaking 


b + sin 6 


', 0< & < 1 


2H-breaking 
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FIG. 2: (Color online) Examples of stable even and odd modes, found in the system with self-focusing nonlinear potential JTJl), 
at rotation speed uj = 0.25, with total power -Peven — -Podd — 

2. Panels (a), (c) display, severally, real and imaginary parts of 
the even and odd modes, while (b), (d) show their local-power (density) profiles. The dashed curves in (b) and (d), and in 
similar panels displayed below, depict the corresponding nonlinearity- modulation profile, cr(6')\ in the present case, it is cos 2 6' . 

III. NUMERICAL RESULTS 

A. The classification of trapped modes 

For all types of the nonlinear potentials denned in Eq. (pQ), the numerical solution of stationary equation ([7]) with 
different inputs (initial guesses) makes it possible to identify five basic species of stationary states, which are listed in 
Table U 

Following the symmetry of the nonlinearity- modulation patterns in Eq. (pQ), the modes are identified as symmetric 
(alias even) and antisymmetric (alias odd), with respect to the central point, 0' = 0. The mode of the "second- 
harmonic" (2H, also even) type refers to the dominant term in its Fourier decomposition. The names of the last two 
modes in Table U come from the types of their symmetry breaking. In particular, the named of 2H-breaking mode is 
come from its spontaneous symmetry breaking via the 2H mode, which will be demonstrated in the following. Note 
that all the inputs displayed in the table are real functions, while the numerically found solutions of Eq. (pp) with 
uo 7^ are actually complex. Naturally, real parts of the solutions generated by the real inputs indicated in the table 
have the same parity (even/odd) as the inputs, while the imaginary parts feature the opposite parity (odd/even) (the 
input of the 2H-breaking type does not feature a certain parity). It is shown below too that maxima and minima of 
the local power (density) of the even mode defined in the table coincide with the maxima and minima of the local 
nonlinearity, cr(#'), while for the odd and 2H modes the relation is the opposite, with the peak powers sitting in 
potential wells [see Eq. (pQ)], hence the odd and 2H modes tend to have lower values of the energy, and may play the 
role of the ground state, as confirmed below. 

B. The self- focusing nonlinearity 

The SF nonlinear potential, defined as per Eq. (pji), gives rise to four types of dynamically stable trapped modes, 
viz., even (symmetric), odd (antisymmetric), 2H and 2H-breaking ones (i.e., only the symmetry-breaking species is 
missing in this case). Typical examples of these stable modes are displayed in Figs. [2] and 02 In addition to the 
above-mentioned fact that the maxima and minima of the local power coincide with those of the nonlinear potential 
for the even mode, and, on the contrary, coincide with minima and maxima of the potential for the odd and 2H modes, 
suggesting that either of the latter modes may be a ground state, the figures demonstrate that the 2H-breaking mode 
has one maximum and one minimum of the local power, both sitting in nonlinear-potential wells (minima). 

Results of the numerical analysis for the stability of the modes in the model with the SF nonlinear potential are 
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(a) (b) (c) (d) 



FIG. 3: (Color online) Examples of stable 2H and 2H-breaking modes, found in the system with the self-focusing nonlinear 
potential at uj — 0.25 and P2H = 1 and P2H-break. = 2, respectively. Panels have the same meaning as in Fig. [2] 




(a) (b) (c) 

FIG. 4: (Color online) (a) and (b): Stability diagrams for the even and odd (antisymmetric) modes, respectively, obtained in 
the system with of the self- focusing nonlinearity, in the plane of the rotation speed (uj) and total power (P). (c) The stability 
diagram for the set of the second-harmonic (2H) and 2H-breaking modes. In panels (a), (b), and (c), respectively, the even, 
odd, and 2H modes are stable in the red areas, and unstable in the blank ones. In panel (c), the 2H-breaking mode is stable in 
the green (middle) area, and the bistability, i.e., coexistence of 2H and 2H-breaking stable modes, occurs in the yellow (right 
edge) region. In the blank area of panel (c), no 2H-breaking mode, stable or unstable one, is found. 



summarized in Fig. 2J in the form of diagrams drawn in the plane of (P,uj). They show that the stability regions 
of the even, odd, 2H and 2H-breaking modes strongly overlap between themselves. In particular, the asymmetry 
parameter of the 2H-breaking modes, which is defined as 



ASP EE 



/ \<t>\ 2 d0'- / 

JO J -it 



IP, (13) 



is displayed in Figs. |5(a)|an d |5(b)| The plots demonstrate that the transition between 2H and 2H-breaking modes is 
of the supercritical type 39j. 

The multistability, which is obvious in Fig. 2J makes it necessary to compare energies of the coexisting dynamically 
stable modes, defined by Eq. ([5]), in order to identify the ground state that realizes the energy minimum. First, in 
Fig. EKa) we show the results along horizontal cuts of all the three panels of Fig. 2J made at a constant value of the 
total power, P = 0.5, with the rotation speed varying in the interval of < uj < 0.4. The tristability of the even, odd 
and 2H modes takes place along this segment. Further, Fig. [6fb) displays the results along vertical cuts of panel (c) 
of Fig. [4] made at uj = 0.5, while the power is varying as 0.2 < P < 1.6 [it is seen in panels (b) and (c) of Fig. [4] 
that the 2H-breaking mode also coexists with the stable odd one, but the energy of the odd mode is definitely larger, 
therefore it is not displayed in Fig. E^b)]. In panel (b) of Fig. [6j the two branches merge into one, with H — » 0, at 
P — )• 0, as the nonlinear potential vanishes in this limit. 

From Fig. 6(a)| we conclude that i^2H < Hodd < H even , while Fig. |6(b)| shows that #2H- break < ^2H- Calculations 



of the energy, performed along other horizontal and vertical cuts, demonstrate that the following relation between the 
energies of the different modes, suggested by these inequalities, is always correct: 

#2H-break < #2H < #odd < H even . (14) 

Thus, the 2H-breaking mode, if it exists [recall that, in the plane of (oj,P) shown in Fig. ffl(b), it exists above the 
curve separating the bottom (red) and middle (yellow) areas], represents the ground state in the system with the 
self- focusing nonlinearity. If the latter mode does not exist, then the 2H state plays the same role. 
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FIG. 5: (Color online) (a)The ASP of the 2H and 2H-breaking modes, defined as per Eq. (|13[) , as a function of the total power, 
P, with the rotation speed fixed at uj = 0.25, in the system with the self-focusing nonlinearity (b). The same as a function of uj, 
at a fixed total power, P = 1. The plot includes regions symmetric with respect to uj — 0.5, to stress the respective symmetry 
of modes in the present system, (c) The ASP of even and symmetry-breaking mode as a function of P for uj = 0.45, in the case 
of the self-defocusing nonlinearity 



0.4- 


— 


0.00- 


0.3 




^.u — 
A— A— A— A— A— A- 


— □ — even 
— O— odd 

P=0.5 


-0.05- 
-0.10 
^ -0.15- 


0.1- 
0.0- 


_A— A— A 


-0.20 
-0.25- 



0.0 0.1 



0.2 
CO 



0.3 0.4 




0.0 0.4 



0.8p 1.2 




0.5 1.0 1.5p2.0 2.5 3.0 



(a) 



(b) 



(c) 



FIG. 6: (Color online) (a) Energies of the even, odd, and second-harmonic modes in the system with the self-focusing nonlin- 
earity, computed, as per Eq. (0, for P = 0.5 and < uj < 0.4. (b) Energies of the 2H-breaking and second-harmonic modes, 
computed for uj = 0.5 and 0.2 < P < 1.6. (c) Energies of the even and odd modes along the vertical segment, with uj = and 
0.5 < P < 3. 



In Fig. [6(a)| we can see that in the limit of uj = 0, the energies of the even and odd modes are very close, which 
seems in contradiction with the fact that these two modes are essentially different, having the opposite parities. The 
reason is that the total power that we chose here (P = 0.5) is not large enough to show the energy difference between 
the two modes. To clarify the point, in Fig. |6(c) we display the energy curves for these two modes selected along the 
vertical segment with uj = and 0.5 < P < 3. It shows that the energy difference indeed increases with the growth 
of P. 



C. The self-defocusing nonlinearity 

In the case of the SDF nonlinear potential, represented by Eq. ([Dd), the numerical solution of Eq. ([7]) reveals stable 
modes of three types, viz., even, symmetry-breaking, and 2H (recall that symmetry-breaking modes were not found in 
the system with the SF nonlinearity), while odd and 2H-breaking states do not exist in this case. Because the profiles 
of the even and 2H modes are quite similar to their counterparts presented above in Fig. [2j we here display, in Fig. [7| 
only a typical stable symmetry-breaking mode. The shape of this mode seems symmetric, centered at 0' = 0; however, 
it is classified as an asymmetric mode, as the power profile of a true symmetric state would be double-humped, cf. 
Figs. E^b,d) and[3^b), while the present one has a single maximum, similar to the intensity distribution in the 2H- 
breaking state in Fig. [3fd). Accordingly, the ASP (effective asymmetry measure) for the symmetry-breaking mode is 



(a) 



(b) 



FIG. 7: (Color online) A stable symmetry-breaking mode, found in the system with the self-defocusing nonlinearity, at (cj, 
(0.45, 1): (a) its real and imaginary parts; (b) the power profile. 
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FIG. 8: (Color online) (a) The stability diagram for the even and symmetry-breaking modes in the system with the self- 
defocusing nonlinearity. The red (largest) and yellow (smallest, near the right edge) areas designate, severally, the stability 
region of the even mode, and the region of the coexistence (bist ability) of the even and symmetry-breaking modes. In the blank 
area, no stable modes of these types are found (in fact, an unstable even mode exists in that area). The second-harmonic mode 
exists and is stable in the entire plane, (b) Energies of the even, symmetry-breaking, and second-harmonic modes along the 
vertical cut of panel (a) at uj — 0.5 and 0.2 < P < 1. 



introduced as follows, instead of the above definition ([T3]) : 
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(15) 



to stress the lack of the asymmetry between the central and peripheral parts parts of the mode. 

The stability and energy diagrams for the even, symmetry-breaking, and 2H modes are displayed in Fig. [8] [the 2H 
mode exists and is stable in the entire plane of (P, a;), therefore it is not specially marked in panel [5Ja)]. In particular, 
it is observed that the symmetry-breaking mode exists near the right edge of the rotational Brillouin zone, i.e., it does 
not exist in the stationary system (with uj = 0). The ASP of the even and symmetry-breaking modes, defined per Eq. 
([T5]h is displayed in Fig. 5(c)| as a function of of the total power. The absence of the linkage between the branches 
representing these two modes implies that they are not related by any bifurcation. 

Figure |8(b)| displays the comparison of energies of these three kinds of the modes (even, symmetry-breaking, and 
2H ones) along the vertical cut made at uj = 0.5, with the power varying in interval 0.2 < P < 1. It demonstrates 
that the curve for the 2H mode goes close to but slightly lower than its counterpart for the symmetry-breaking mode. 
The analysis of more general data, produced by the numerical calculations for the system with the SDF nonlinear 
potential, demonstrates that energies of all the three dynamically stable modes existing in this case are ordered as 
follows: 



i^2H < ^symm- break < H ( 



even? 



(16) 



cf. Eq. (Q3I). Thus, the 2H mode plays the role of the ground state in the case of the SDF nonlinearity (recall this 
mode exists at all values of uj and P). 
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(a) (b) (c) 



FIG. 9: (Color online) Stability diagrams for the system with the alternating self-focusing - self-defocusing nonlinearity. (a) 
Even and symmetry-breaking modes; (b) odd and 2H-breaking modes; (c) the 2H modes. In (a), the even mode is stable in the 
red (bottom) area; in the small yellow region it coexists with the symmetry-breaking one, and no stable modes of these types 
are found in the blank area. In (b), the odd mode is stable in the red (bottom) area, the 2H-breaking mode is stable in the 
green region (adjacent to the right edge of the panel), and these two types coexist in the yellow area. In (c), the 2H mode is 
stable in the red area. 

D. The alternating self-focusing — self-defocusing nonlinear potential 

In the case of the alternating SF-SDF nonlinearity, defined as per Eq. (pQ), the numerical solutions reveal the 
existence of all the five types of stable trapped modes indicated in Table [U Profiles of these modes are quite similar 
to those of their counterparts displayed above in Figs. [2j [3j and [7J therefore we do not show them again here. The 
respective stability diagrams in the (P, lj) plane are presented in Fig. [9] In particular, the absence of a bistability 
area in panel (c) of Fig. [9] suggests that the transition between the 2H and 2H-breaking modes is supercritical, like 
in the system with the SF nonlinearity, cf. Fig. ??. 

The comparison of the energies of the five types of the modes which may be stable in the case of the sign-alternating 
nonlinearity demonstrates the following ordering, cf. Eqs. JT2J) and ([T6]) : 

#2H- break < #2H < H sy mm. -break, #odd < H even . (17) 

The energies of the symmetry-breaking and odd modes are not compared in Eq. (fT7|) . as their stability regions do 
not overlap, see Fig. [9]^a,b). Thus, the 2H-breaking mode, when it exists, plays the role of the ground state in the 
present case; otherwise, the ground state is represented by the 2H mode, see panels (b) and (c) in Fig. [9j 

IV. THE ANALYTICAL APPROACH 

A. The two-mode approximation 

The present setting may be naturally approximated by a finite-mode truncation of the expansion of stationary field 
4>{x) over the set of spatial harmonics. The simplest approximation reduces to the substitution of truncation 

0(0') = a + ai exp (iO') (18) 

into Eq. ([7]). In this expression, ao may be fixed to be real, while amplitude a\ is allowed to be complex. This 
approach is consistent if, in the linear approximation, each term in combination ([T8]) is an exact solution of Eq. ([7]) 
for a common value of /i, hence the zeroth-order approximation exists, and the analysis of weakly nonlinear states 
can be developed around it. It is easy to see that such a case corresponds, in the zeroth approximation, to uj = 1/2 
[which is exactly the right edge of zone (fit)]) ] and \i = [34[. Then, weakly nonlinear modes can be constructed in an 
approximate analytical form, assuming that /i, a^, |«i| 2 , and 

5=l/2-a; (19) 

are all small quantities. To this end, ansatz (fT8j) and a particular expression (pQ) for a (9 r ) are substituted into Eq. 
0, and equations for amplitudes ao and a\ are derived as balance conditions for the zeroth and first harmonics. 
Ansatz (fT8|) corresponds to the following approximation for the total power (j4]), 



P= l + |</>|W = 2^(ag + |a 1 | 2 ), 

J — 7T 



(20) 
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which will be used below too. 



B. The self- focusing nonlinearity 

In the case of the SF nonlinear potential (pji), the two- mode approximation ([T8]) leads to the following equations: 

A* = ~2 a o - M + 4 a i> 

M _ 5= _ a 2_l| ai |2 + 1 2^1 (21) 
Z 4 (2i 

As shown above by the numerical analysis, the SF nonlinearity gives rise, inter alia, to the 2H-breaking mode, which 
is generated from input b + sin^ in Table HI To capture the part of the solution corresponding to sin^ in the input, 

we set 

a x = ic, (22) 

Then, Eq. Q2D yields 

ag = (4/21) (5* -3/x), 

c 2 = -(4/21) (2(5 + 3/i). (23) 

Substituting solution ([23]) into expression ([20]) . we obtain a relation between the propagation constant and total power 
for the 2H-breaking mode, 

fi = - - —P. (24) 
P 2 16tt V ; 

It demonstrates that this mode, as predicted by the analytical approximation, satisfied the Vakhitov-Kolokolov (VK) 
criterion, dfi/dP < 0, which is a necessary stability condition for patterns supported by the SF nonlinearity [40|. 
Furthermore, Eq. ([23]) predicts the emergence of the 2H-breaking modes at Oq = 0, i.e., 

fi = 55/3. (25) 

Since the modes with rotation speeds related by Eq. (j9]) are mutually tantamount, Eqs. ([25]) and ([24]) predict the 
coexistence of the odd and 2H-breaking modes at 

P>(Pmm) = (8ir/3)\5\. (26) 

This analytical result is compared with the corresponding numerical findings in Fig. fTUT a). which shows a reasonably 
good agreement. 



C. The self-defocusing nonlinearity 

The substitution of the same ansatz ([18]) into Eq. ([7j) in the case of the SDF nonlinear potential, taken as per Eq. 
(jib), yields the following algebraic equations, instead of Eq. ([2T]) derived above for the SF nonlinearity: 

1 o I i o 1 2 

fi= -o +|ai| +^%, 

f 1 -S = a 2 + ha 1 \ 2 + \a 2 ^-. (27) 
z 4 a\ 

The above numerical results for the SDF nonlinearity demonstrate the existence of the symmetry-breaking mode in 
this case, which is generated by input 6 + cos^ in Table [U To capture this mode by means of ansatz JT5J), it is natural 
to set 



ai = c = real, 



(28) 



10 



on the contrary to Eq. ([22]) , where a\ was imaginary. In this case, the algebraic system ([27]) yields 

al = (4/21) (3/X-55), 

c 2 = (4/21)(3/x + 2(J), (29) 

cf. Eq. ([23]) . According to Eq. ([20]) . the relations between the propagation constant and total power takes the 
following form for solution ([29]) : 

cf. Eq. ([M]) . This equation demonstrates that the asymmetric mode satisfies the anti-VK criterion, dfi/dP > 0, 
which, as argued in Ref. |41[, may play the role of a necessary stability condition for modes supported by the SDF 
nonlinearity. 

The analytical approximation predicts the existence boundary for the symmetry-breaking states, al = 0, in the 
same form ([25]) as it was obtained above for the SF nonlinearity. Consequently, the existence region for these states is 
predicted in the form coinciding with that given by Eq. ([26]) . This result is compared with its numerical counterpart 
in Fig. fTm». 



The corresponding relation between \i and P is 



D. The alternating self-focusing — self-defocusing nonlinearity 

In the case of the alternating SF-SDF nonlinear potential ([It), the substitution of ansatz ([T8]) into Eq. ([7]) yields 
the algebraic equations in the form which is somewhat simpler than Eqs. ([2T]) and ([27]) derived above for the "pure" 
SF and SDF nonlinear it ies: 

1 2 

A* = 2 a i> 

H-5 = \alf. (31) 

First, if we choose real ai, as in Eq. ([28]) . which refers to the symmetry-breaking mode, Eq. ([3T]) yields 

al = 2(/i - J), 

c 2 = 2/i. (32) 

which satisfies the above-mentioned anti-VK criterion, cf. Eq. ([30]) . i.e., it is possible to assume that the stability of 
the symmetry-breaking mode is supported by the SDF part of the alternating nonlinear potential. Further, Eq. ([32]) 
predicts the emergence of the symmetry-breaking mode (a 2 , = 0) at /i = 5. In terms of the total power related to /i 
by Eq. ([20]) . this implies that this mode exists at 

P>Pmin = ^\S\, (34) 

cf. Eq. ([26]) . The comparison of this analytical result with its numerical counterpart in shown in Fig. [TUt c). 

The numerical results presented in the previous section demonstrate that the alternating SF-SDF nonlinear potential 
support the stable 2H-breaking mode too. To describe it in the framework of the two-mode approximation, we now 
assume a\ to be imaginary, as in Eq. ([22]) . In this case, Eq. ([32]) yields 

al = 2(6 - /i), 

c 2 = -2/i, (35) 
cf. Eq. ([3T]) . the respective relation between \i and P being 

"-5-S- (36) 

cf. Eq. ([33]) . The latter relation satisfies the VK criterion, which implies that the stability of the 2H-breaking mode is 
supported by the SF part of the alternating nonlinear potential. The emergence of the 2H-breaking mode corresponds 
to al = 0, i.e., again /i = as in the case of the symmetry-breaking mode, under the same alternating nonlinear 
potential. Finally, this means that the existence of the stable 2H-breaking mode is predicted in the same region ([34]) 
as for its symmetry-breaking counterpart. The latter prediction is compared to the numerical findings in Fig. [TOT d). 
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FIG. 10: (Color online) The comparison between the theoretically predicted (dashed lines) and numerically found boundaries 
of the existence of the stable 2H-breaking (a,d) and symmetry-breaking (b,c) modes in the plane of (S,P), in the interval of 
< \8\ = 1 1/2 — uj\ < 0.05. In the red areas, only even or 2H modes are produced by the numerical solution, while in the yellow 
(bistability) regions they coexist with the stable symmetry- or 2H-breaking modes, respectively. Panels (a) and (b) pertain, 
severally, to the SF and SDF nonlinear potentials, while (c) and (d) correspond to the alternating SF-SDF potential. These 
four panels are, actually, zoomed versions of the right bottom corners of Figs. |4(b)| |8(a)| |9(a)| and |9(b)| respectively, with 
analytically predicted lines ([26]) or (|34|) added to each panel. 



V. CONCLUSION 

This work aimed to study the existence and stability of modes trapped in the rotating nonlinear-lattice potentials 
of the SF, SDF (self-focusing and defocusing) and alternating SF-SDF types. The consideration was carried out for 
the first rotational Brillouin zone in the rotating reference frame, i.e., for < uj < 1/2, where uj is the rotation speed. 
The stability analysis was performed through the computation of eigenvalues for small perturbations, and verified by 
direct simulations. The model can be realized in the spatial domain, in terms of a twisted-pipe optical waveguide, with 
the built-in azimuthal modulation of the local Kerr coefficient, or, in the temporal domain, as the Gross-Pitaevskii 
equation for BEC loaded into a toroidal trap, under the action of a rotating optical or magnetic structure which 
affects the local value of the scattering length. 

In the SF system, four types of different modes have been identified: even, odd, which are dominated by combinations 
of the fundamental and zeroth angular harmonic, and the 2H (second-harmonic) and 2H-breaking states. On the other 
hand, the SDF nonlinear potential supports three species of the trapped states: 2H, even, and symmetry-breaking 
ones, the latter existing only at the rotation speed close to the right edge of the Brillouin zone, uj = 1/2, and in a 
limited interval of values of the total power, P. As for the alternating SF-SDF nonlinear potential, it supports all 
the five species of the trapped modes. Transitions between the 2H and 2H-breaking modes are of the supercritical 
types. The energy comparison reveals that, in the SF and SF-SDF systems alike, the 2H-breaking mode, if it exists, 
represents the ground state; otherwise, this role is played by the 2H mode. The ground state of the SDF system is 
always represented by the 2H solution. 

This work may be naturally developed in other directions. On the one hand, it is relevant to consider rotating 
nonlinear lattices with smaller azimuthal periods, 27r/n, for integer n > 1, unlike the case of n = 1 investigated here. 
On the other hand, it may be interesting to consider a two-dimensional version of the present one-dimensional model 
(with an entire rotating plane, rather than the thin ring, which will include effects of the Coriliolis' force). 
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